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Abstract: 

The  long  term  behaviour  of  tropical  cyclones  in  the  prototype  problem  for  cyclone  intensification  on  an  /-plane  is  examined  using  a 
nonhydrostatic,  three-dimensional,  numerical  model.  After  reaching  a  mature  intensity,  the  model  storms  progressively  decay  while 
both  the  inner  core  size,  characterized  by  the  radius  of  the  eyewall,  and  size  of  the  outer  circulation,  measured  for  example  by  the  radius 
of  gale  force  winds,  progressively  increase.  This  behaviour  is  explained  in  terms  of  a  boundary-layer  control  mechanism  in  which  the 
expansion  of  the  swirling  wind  in  the  lower  troposphere  leads  through  boundary-layer  dynamics  to  an  increase  in  the  radii  of  forced 
eyewall  ascent  as  well  as  to  a  reduction  in  the  maximum  tangential  wind  speed  in  the  layer.  These  changes  are  accompanied  by  changes 
in  the  radial  and  vertical  distribution  of  diabatic  heating.  As  long  as  the  aggregate  effects  of  inner-core  convection,  characterized  by 
the  distribution  of  diabatic  heating  are  able  to  draw  absolute  angular  momentum  surfaces  inwards,  the  outer  circulation  will  continue 
to  expand.  The  quantitative  effects  of  latitude  on  the  foregoing  processes  are  investigated  also. 

The  study  provides  new  insight  on  the  factors  controlling  the  evolution  of  size  and  intensity  of  a  tropical  cyclone.  It  provides  also  a 
plausible,  and  arguably  simpler,  explanation  for  the  expansion  of  the  inner-core  of  Hurricane  Isabel  (2003)  than  that  given  previously. 
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1  Introduction 

When  measured  by  the  extent  of  gale-force  winds  (wind 
speeds  exceeding  17  m  s-1),  tropical  cyclones  are  observed 
to  possess  a  wide  range  of  sizes  from  large  Pacific  typhoons 
with  gale  force  winds  extending  beyond  1000  km  from  their 
centre  to  midget  storms  with  a  radius  of  gales  no  more  than 
about  50  km.  There  is  a  large  variation  also  in  the  size 
of  the  eye,  which  is  not  necessarily  related  to  the  extent 
of  gales.  Examples  of  such  extreme  sizes  include  at  the 
large  end  supertyphoon  Tip  (1979),  which  occurred  in  the 
western  North  Pacific  (Dunnavan  and  Diercks,  1980),  and 
hurricane  Sandy  (2012)  which  caused  widespread  damage 
along  the  eastern  seaboard  of  the  United  States  (e.g.  Lussier 
et  ah,  2015).  An  example  at  the  small  end  is  tropical 
cyclone  Tracy  (1974),  which  devastated  the  Australian  city 
of  Darwin  on  Christmas  Day  (Davidson  2010  and  refs).  The 
radius  of  gales  in  the  case  of  Tracy  was  merely  50  km,  less 
than  the  size  of  the  eye  in  Tip,  where  gales  were  spread  over 
a  1 100  km  radius.  Thus  Tip  had  an  area  about  400  hundred 
times  larger  than  that  of  Tracy.  In  contrast,  hurricane  Wilma 
(2005)  was  a  case  with  a  broad  tangential  circulation,  but  a 
pinhole  sized  eye  (Beven  et  ah,  2008). 

In  a  climatological  study  of  tropical  cyclone  size,  Mer¬ 
rill  (1984)  showed  that  tropical  cyclones  of  the  western 

1  Correspondence  to:  Dr.  Gerard  Kilroy.  Meteorological  Institute, 
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North  Pacific  are  characteristically  twice  as  large  as  their 
Atlantic  counterparts,  and  that  the  typical  size  of  tropi¬ 
cal  cyclones  varies  seasonally  and  regionally  and  is  only 
weakly  correlated  with  cyclone  intensity  (as  measured  by 
the  maximum  surface  wind  or  minimum  surface  pressure). 
He  noted  that  within  a  single  season,  tropical  cyclones 
come  in  all  sizes  and  attain  a  variety  of  intensities.  A  more 
recent  climatological  study  of  tropical  cyclone  size  by  Dean 
et  ah  (2009)  showed  that  the  size  distribution  of  Atlantic 
tropical  cyclones  was  log-normal1. 

Although  there  has  been  a  number  of  theoretical 
studies  examining  factors  that  determine  tropical  cyclone 
size  (e.g.  Yamasaki,  1968;  Rotunno  and  Emanuel,  1987; 
DeMaria  and  Pickle,  1988;  Smith  et  ah,  2011;  Rappin  et  ah, 
201 1 ;  Li  et  ah,  2012;  Hakim,  2011;  Chavas  and  Emanuel, 
2014),  there  remain  some  basic  issues  to  be  resolved.  For 
example,  an  underlying  assumption  of  all  of  these  studies  is 
that  there  exists  a  global  quasi-steady  solution  for  storms, 
an  assumption  that  would  require,  inter  alia,  that  the  storm 
environment  be  quasi-steady.  However,  the  existence  of  a 
realistic  globally,  quasi-steady  state,  even  when  the  envi¬ 
ronment  is  quiescent,  has  been  questioned  on  the  grounds 

hThe  metric  for  size  used  by  these  authors  was  “the  radius  of  vanishing 
storm  winds  normalized  by  the  theoretical  upper  bound  for  this  radius 
given  by  the  ratio  of  the  potential  intensity  to  the  Coriolis  parameter”.  The 
radius  of  vanishing  storm  winds  was  inferred  from  the  estimated  radius  of 
gale  force  winds  using  Emanuel’s  revised  steady-state  hurricane  model 
(Emanuel,  2004). 
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of  both  angular  momentum  and  thermodynamic  considera¬ 
tions  (Smith  et  al.  2014),  raising  the  likelihood  that  both 
the  size  and  intensity  of  the  storm  evolve  progressively 
with  time.  In  fact,  according  to  the  conventional  paradigm 
for  tropical  cyclone  intensification  (see  Montgomery  and 
Smith,  2014,  and  references),  one  would  anticipate  that 
the  outer  circulation  will  expand  as  long  as  the  aggregate 
effect  of  deep  convection  (including  the  eyewall  and  con¬ 
vective  rainbands  -  Fudeyasu  and  Wang,  2011)  remains 
strong  enough  to  maintain  the  inward  migration  of  abso¬ 
lute  angular  momentum  surfaces.  If  this  is  the  case,  the 
broadening  circulation  would  have  consequences  for  the 
boundary  layer  dynamics,  which  play  a  role  in  determining 
the  radii  at  which  air  ascends  into  the  eyewall  and  the  maxi¬ 
mum  tangential  wind  speed,  which  has  been  shown  to  occur 
in  the  boundary  layer  (e.g.  Kepert  2006a,b;  Montgomery 
et  al.  2006,  2014;  Bell  and  Montgomery  2008;  Sanger  et  al. 
2014). 

Xu  and  Wang  (2010)  suggested  that  a  broadening 
of  the  circulation  would  increase  also  the  surface  mois¬ 
ture  flux  outside  the  eyewall.  thereby  increasing  the  “con¬ 
vective  activity”  at  these  radii  and  enhancing  the  radial 
inflow  beyond  these  radii.  While  at  first  sight  this  argu¬ 
ment  may  seem  reasonable,  it  is  important  to  quantify  the 
link  between  the  increase  in  “convective  activity”  as  a  result 
of  moistening.  Presumably,  “convective  activity”  is  charac¬ 
terized  by  the  convective  mass  flux  and  its  variation  with 
radius. 

In  two  recent  papers,  Chan  and  Chan  (2014)  and 
Frisius  (2015)  carried  out  a  series  of  idealized  numerical 
simulations  of  tropical  cyclones  with  a  view  to  exploring 
the  controls  on  storm  size.  Chan  and  Chan  (2014)  sought 
“  ...  to  understand:  (our  emphasis)  how  the  initial  size  ... 
and  planetary  vorticity  ...  influence  TC  (tropical  cyclone: 
our  insertion)  size  change”.  They  defined  the  size  as  the 
radius  of  the  azimuthally-averaged  wind  speed  of  17  m 
s-1  at  a  height  of  10  m.  An  understanding  of  size  using 
this  definition  requires  a  consideration  of  frictional  effects, 
but  no  attention  was  provided  in  this  regard.  Instead,  size 
differences  were  explained  mostly  in  terms  of  differences 
in  inertial  stability,  the  relevance  of  which  is  questionable 
in  the  friction  layer2.  In  addition,  they  did  not  mention  the 
differences  in  the  strength  of  the  convective  forcing,  which, 
as  we  will  show  here,  are  fundamental  to  understanding 
these  effects.  Indeed,  Smith  et  al.  (2015)  showed  that  the 
convective  forcing  is  paramount  and  that  inertial  stability 
(above  the  boundary  layer)  plays  a  minor  role,  certainly  on 
the  intensification  at  different  latitudes.  They  showed  also 
that  the  boundary  layer  has  a  strong  control  on  the  location 
and  strength  of  the  convective  forcing  as  the  latitude  is 
varied. 

2The  point  is  that  inertial  stability  is  normally  judged  by  considering  the 
forces  on  an  air  parcel  that  is  subjected  to  a  small  radial  displacement  in 
a  swirling  flow  in  gradient  wind  balance.  Such  an  analysis  is  meaningless 
in  the  boundary  layer,  where  gradient  wind  balance  is  not  satisfied  and 
where,  at  outer  radii,  the  agradient  force  is  negative. 


A  significant  finding  of  Frisius  op.  cit.,  who  used 
an  axisymmetric  model,  was  that,  even  in  a  quiescent 
environment,  a  steady  state  solution  of  the  model  requires 
an  artificial  (and  therefore  unrealistic)  source  of  absolute 
angular  momentum  to  maintain  the  angular  momentum 
lost  by  the  frictional  torque  near  the  surface  where  the 
flow  is  cyclonic.  When  this  source  was  suppressed,  the 
upper  anticyclone  descended  to  the  surface,  where  cyclonic 
relative  angular  momentum  can  be  diffused  into  the  system. 
Both  these  results  are  in  accordance  with  the  theoretical 
predictions  of  Smith  et  al.  (2014).  Frisius  (2015,  pl4) 
conceded  that  his  steady  state  solutions,  which  took  many 
months  to  achieve,  were  not  relevant  to  understanding  size 
in  realistic  storms,  which  are  intrinsically  unsteady. 

Chavas  and  Emanuel  (2014)  used  extended  time  inte¬ 
grations  of  an  axisymmetric  model  to  examine  the  size  of 
tropical  cyclones  in  a  state  of  statistical  equilibrium.  They 
found  that  the  size  scales  linearly  with  potential  intensity 
divided  by  the  Coriolis  parameter.  Like  Frisius,  they  con¬ 
ceded  (pi 678)  that  “the  extent  to  which  these  equilibrium 
results  can  be  applied  to  real  storms  in  nature  is  not  clear”, 
and  they  explain  the  reasons  why:  “the  time  scales  to  equi¬ 
librium  identified  here  for  the  control  simulation  are  sig¬ 
nificantly  longer  than  the  lifespan  of  tropical  cyclones  on 
Earth”,  and  “storms  in  nature  rarely  exist  in  a  truly  quasi¬ 
steady  environment  for  more  than  a  couple  of  days,  if  at 
all.” 

We  acknowledge  the  challenges  of  understanding  size 
change,  but  do  not  subscribe  to  Frisius’s  statement  on  the 
first  page  of  his  paper  that  “little  is  known  about  the  mech¬ 
anism  controlling  (tropical  cyclone:  our  insertion)  size”. 
Nevertheless,  we  think  that  there  is  scope  for  improving  an 
understanding  of  the  interplay  of  processes  that  must  be  of 
prime  importance  and  such  is  the  motivation  of  the  present 
paper.  Any  improvement  would  have  obvious  benefits  for 
forecasters,  if  only  in  improving  their  ability  to  interpret 
model  behaviour,  on  which  all  modern  forecasts  are  based. 

The  present  paper  builds  on  the  studies  of  Smith  et  al. 
(2014)  and  Smith  et  al.  (2015,  -  henceforth  SKM).  These 
focused  on  the  prototype  problem  for  tropical  cyclone 
intensification,  which  considers  the  evolution  of  a  pre¬ 
scribed,  initially  cloud-free  axisymmetric  vortex  in  a  qui¬ 
escent  environment  on  an  /-plane.  The  first  of  these  studies 
sought  to  understand  the  lack  of  a  steady  state  in  the  30  day 
simulation,  while  the  second  sought  to  understand  why  low 
latitude  storms  intensify  more  rapidly  that  those  at  higher 
latitudes,  all  other  factors  being  the  same.  Here  we  use  the 
same  three-dimensional  numerical  model  as  SKM  to  under¬ 
stand  the  long  term  behaviour  of  the  vortex  (on  the  time 
scale  of  a  month)  and  the  effects  of  latitude  on  this  evolu¬ 
tion. 

The  paper  is  organized  as  follows.  In  section  2  we 
present  the  configuration  of  the  numerical  model,  then 
in  section  3  we  discuss  the  results  of  the  calculations. 
In  section  4  we  provide  physical  interpretations  of  the 
model  behaviour  and  in  section  5  we  offer  an  alternative 
explanation  for  size  expansion  in  two  notable  tropical 
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cyclones.  In  section  6  we  discuss  the  latitudinal  dependence 
of  the  calculations.  The  conclusions  are  given  in  section  7. 


2  The  numerical  model 

The  numerical  experiments  use  a  modified  version  of  the 
Pennsylvania  State  University/National  Center  for  Atmos¬ 
pheric  Research  mesoscale  model  (MM5)  and  are  with  a 
few  exceptions  identical  to  those  described  by  SKM:  see 
their  section  2.  A  detailed  description  of  the  MM5  model 
can  be  found  in  Grell  et  al.  (1995).  The  model  is  config¬ 
ured  here  with  three  domains:  a  coarse  mesh  of  45 -km  grid 
spacing  and  two,  two-way  nested  domains  of  15  and  5  km 
grid  spacing,  respectively.  The  domains  are  square  and  are 
9000  km,  4500  km,  1500  km  on  each  side.  The  model  has 
24  er-levels  in  the  vertical,  7  of  which  are  below  850  mb  to 
provide  adequate  vertical  resolution  of  the  boundary  layer. 
A  key  difference  between  the  model  configuration  here  and 
that  used  by  SKM  is  the  application  of  Newtonian  relax¬ 
ation  to  the  temperature  field  with  a  timescale  of  10  days 
(Mapes  and  Zuidema  1996). 

The  initial  vortex  is  axisymmetric  with  a  maximum 
tangential  wind  speed  of  15  m  s_1  at  the  surface  at  a  radius 
of  100  km.  The  magnitude  of  the  tangential  wind  decreases 
sinusoidally  with  height,  vanishing  at  the  top  model  level. 
The  temperature  field  is  initialised  to  be  in  thermal  wind 
balance  with  the  wind  field  using  the  method  described  by 
Smith  (2006).  The  far-field  temperature  and  humidity  are 
based  on  the  neutral  sounding  of  Rotunno  and  Emanuel 
(1987). 

Three  calculations  are  performed  on  an  /-plane  cen¬ 
tred  at  latitudes  10°N,  20°N  and  30°N.  The  simulations  are 
carried  out  for  30  days  with  data  output  every  15  min.  As 
the  simulations  run  for  30  days,  the  latent  and  sensible  heat 
fluxes  are  switched  off  beyond  a  radius  of  500  km  as  a  pre¬ 
caution  to  avoid  excessive  low-level  moistening  of  the  far 
field  environment.  This  is  another  difference  to  the  calcula¬ 
tions  in  SKM. 


3  Results:  vortex  evolution 

The  early  evolution  of  the  model  cyclone  is  the  same  as  that 
described  in  Nguyen  et  al.  (2008).  In  brief,  the  imposition 
of  surface  friction  from  the  initial  instant  leads  to  the 
formation  of  a  surface-based  boundary  layer  with  inflow 
near  the  surface,  ascent  inside  the  radius  of  maximum 
tangential  wind  speed,  and  outflow  in  a  shallow  layer  above 
the  inflow  layer.  The  inflow  progressively  moistens  through 
surface  moisture  fluxes  and  after  about  10  h,  some  of  the 
ascending  air  reaches  its  level  of  free  convection  and  cells 
of  deep  convection  develop.  At  first,  a  regular  pattern  of 
convective  updraughts  form  in  an  annular  region  inside  the 
initial  radius  of  the  maximum  tangential  wind  speed,  but 
the  pattern  progressively  loses  any  degree  of  symmetry 
reflecting  the  stochastic  nature  of  deep  convection  and 
the  number  of  updraughts  declines.  The  updraughts  rotate 


(a)  time  (d) 


(b)  time  (h) 


(c)  time  (d) 

Figure  1 .  Time  series  of  azimuthally-averaged  quantities  at  latitudes 
10°N,  20°N,  and  30°N:  (a)  maximum  tangential  wind  speed  (Vmax) 
during  the  entire  30  day  calculation  and  (b)  the  first  3  days  which 
encompasses  the  intensification  phase.  Panel  (c)  shows  the  radius  at 
which  the  maximum  tangential  wind  speed  occurs  {Rvm ax)- 


cyclonically  around  the  vortex  centre  and  have  lifetimes 
on  the  order  of  an  hour.  As  they  develop,  they  tilt  and 
stretch  the  local  vorticity  field  and  an  approximate  ring-like 
structure  of  intense,  small-scale,  vorticity  dipoles  emerges. 
The  dipoles  are  highly  asymmetric  in  strength  with  strong 
cyclonic  vorticity  anomalies  and  much  weaker  anticyclonic 
vorticity  anomalies. 

The  development  of  the  updraughts  heralds  a  period 
lasting  about  three  to  six  days  during  which  the  vortex 
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intensifies  rapidly.  Eventually,  the  maximum  azimuthally- 
averaged  tangential  wind  speed,  Vmax,  reaches  a  quasi¬ 
steady  state.  However,  in  terms  of  other  metrics,  this  quasi¬ 
steady  state  is  an  illusion  and  other  quantities  such  as  the 
strength  of  the  upper  anticyclone  and  the  radius  of  gales 
are  still  evolving  (see  e.g.  Smith  et  al.  2014). 

The  longer  term  behaviour  of  the  vortex  is  summarized 
in  Figs.  1  and  2. 

3. 1  Intensity  and  inner-core  size 

Figure  1  shows  time  series  of  Vmax  over  long  and  short 
time  periods  together  with  the  radius  Rvmax  at  which 
Vmax  occurs  in  the  three  calculations  for  different  latitudes. 
These  quantities  can  be  at  any  height,  but  are  typically  at 
a  few  hundred  metres  above  the  surface,  within  the  inflow 
layer.  In  terms  of  Vmax,  the  vortex  in  all  three  simulations 
undergoes  an  intensification  phase  lasting  about  three  to 
six  days,  followed  by  a  quasi-steady  state  which  lasts  up 
to  about  six  days.  All  three  vortices  have  a  similar  Vmax 
during  the  mature  phase.  From  day  nine  onwards,  the 
vortices  progressively  decay,  and  the  decay  rate  is  weaker  in 
the  10°N  simulation  (Fig  la).  At  the  end  of  the  simulation, 
Vmax  in  the  10°N  experiment  has  decreased  by  about  20 
m  s-1  from  its  maximum  value,  while  that  in  the  30°N 
experiment  has  decreased  by  about  35  m  s_1. 

Initially  Rvmax  is  located  at  a  radius  of  100  km 
but  contracts  to  about  20-30  km  after  three  days.  The 
most  rapid  contraction  occurs  during  the  rapid  intensifi¬ 
cation  phase  where  absolute  angular  momentum  surfaces 
are  drawn  inwards  quickly  within,  and  above,  the  bound¬ 
ary  layer.  During  the  most  intense  phase  (between  about 
three  and  nine  days),  Rvmax  is  essentially  the  same  in  all 
three  simulations,  but  it  fluctuates  more  at  20°N  during  the 
decay  phase  (Fig  lc). 

3.2  Outer-core  size 

The  outer  core  size  is  characterized  by  the  outermost 
radius  of  gale  force  winds,  Rgaies ,  defined  here  as  the 
radius  of  17  m  s-1  tangential  winds  at  a  height  of  1 
km,  which  is  approximately  at  the  top  of  the  boundary 
layer.  Forecasters  are  interested  in  the  radius  of  gales,  say 
RgaiesF,  defined  as  the  (outer)  radius  at  which  the  total 
wind  speed  at  a  height  of  10  m  is  17  m  s-1.  However,  from 
a  theoretical  viewpoint,  a  definition  based  on  the  tangential 
wind  component  just  above  the  top  of  the  boundary  layer 
is  a  preferable  measure  of  size,  as  this  radius  is  related 
to  absolute  angular  momentum  at  Rgaies »  ic.  Mgaies  - 
17  x  Rgaies  +  \f  Rgaies ■  Moreover,  above  the  boundary 
layer,  where  frictional  effects  can  be  neglected,  absolute 
angular  momentum  is  materially  conserved  to  a  good 
approximation,  i.e.  dM/dt  =  — us  •  VM,  where  us  is  the 
velocity  vector  in  a  vertical  plane  that  characterizes  the 
secondary  circulation.  In  other  words,  a  spin  up  of  the 
tangential  wind  requires  a  component  of  flow  normal  to  the 
M-surfaces  towards  lower  M. 


For  an  inertially-stable  warm-cored  vortex,  dM/dr  > 
0  and  dM/dz  <  0  so  that  VM  points  downwards  and  out¬ 
wards.  It  follows  that  the  radius  of  gales  above  the  boundary 
layer  will  always  expand  if  the  foregoing  condition  is  satis¬ 
fied.  Clearly,  the  magnitude  of  us  will  be  determined,  inter 
alia,  by  some  measure  of  the  strength  and  radial  distribu¬ 
tion  of  the  convective  forcing,  which,  in  turn,  will  depend 
on  both  the  boundary-layer  dynamics  and  thermodynamics. 
It  may  depend  also  on  the  inertial  stability  above  the  bound¬ 
ary  layer. 

Within  the  boundary  layer,  the  absolute  angular 
momentum  of  an  air  parcel  is  reduced  by  friction.  Thus 
friction  needs  to  be  considered  when  calculating  the  radial 
movement  of  the  M-surfaces  to  determine  whether  the 
radius  of  gales  increases  or  not.  It  follows  that  if  one 
uses  RgaiesF  to  measure  size,  one  cannot  ignore  frictional 
effects  in  any  explanation  of  size  changes.  Moreover,  since 
RgaiesF  is  based  on  the  total  wind  speed  at  a  height  of  10 
m,  the  radial  momentum  equation  must  be  included  in  any 
analysis  of  its  behaviour. 

As  the  latitude  increases,  so  does  the  radial  gradient 
of  absolute  angular  momentum,  but  the  implied  increase 
in  the  inertial  stability  and  any  decrease  in  the  strength  of 
convective  forcing  (see  below),  may  reduce  the  strength  of 
the  secondary  circulation.  Therefore,  the  net  effect  on  the 
radial  displacement  of  the  M-surfaces  and  hence  on  the 
spin  up  rate  cannot  be  anticipated  and  must  be  calculated. 
Figure  2  shows  time  series  of  RgaieS,  RgaiesF  and  the 
difference  between  the  total  winds  at  the  surface  (17  m 
s-1)  at  RgaiesF  and  the  tangential  wind  component  at  a 
height  of  1  km.  As  shown  in  Fig  2a  and  b,  the  outer  core 
size  of  the  vortices  at  10°N  and  20°N  grows  continuously, 
the  most  rapid  growth  occurring  during  the  intensification 
phase.  At  30°N  the  vortex  size  remains  relatively  constant 
during  the  decay  phase  until  the  final  days  of  the  simulation, 
suggesting  that  deep  convection  in  the  weakening  system 
becomes  progressively  unable  to  draw  in  M-surfaces  in  the 
far  field.  Despite  the  larger  background  rotation  at  higher 
latitudes,  the  30°N  vortex  is  the  smallest  in  size  from  day 
15  onwards.  The  20°N  vortex  becomes  the  largest  in  size 
after  about  14  days  and  it  continues  to  grow  in  size  until  the 
end  of  the  simulation,  suggesting  that  — us  •  VM  remains 
positive  at  Rgaies ,  even  though  the  inner-core  intensity  of 
the  vortex  is  decreasing.  In  the  final  two  days  the  10°N 
vortex  becomes  the  largest  in  size,  despite  spending  the 
initial  15  days  as  the  smallest. 

A  comparison  of  the  top  two  panels  in  Fig.  2  shows 
that  the  evolution  of  storm  size  based  on  RgaiesF  is  similar 
to  that  based  on  RgaieS,  although  Rgaies  always  exceeds  the 
value  of  RgaiesF  at  all  three  latitudes.  Typically,  Rgaies  is 
about  25  km  larger  than  RgaiesF  during  the  mature  phase, 
and  the  difference  increases  to  values  as  large  as  60  km  in 
the  decay  phase  at  all  three  latitudes.  Moreover,  at  RgaiesF, 
the  tangential  wind  speed  at  a  height  of  1  km  exceeds  the 
total  surface  wind  speed  by  between  3  and  5  ms-1  with  a 
noticeable  dependence  on  latitude  (Fig.  2c). 
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(a)  time  (d) 


(b)  time  (d) 


(c)  time  (d) 

Figure  2.  Time  series  of  azimuthally-averaged  quantities  at  latitudes 
10°N,  20°N,  and  30°N:  (a)  the  radius  at  which  gale  force  winds 
occurs  ( Rgaies ),  where  Rgaies  is  calculated  at  a  height  of  1  km,  and 
corresponds  to  the  radius  of  17  m  s_1  tangential  winds  outside  the 
eyewall;  (b)  the  radius  at  which  gale  force  winds  occurs  ( RgaiesF ), 
where  RgaiesF  is  calculated  at  a  height  of  10  m,  and  corresponds 
to  the  radius  of  17  m  s-1  total  winds  outside  the  eyewall;  (c)  the 
difference  between  the  tangential  wind  at  a  height  of  1  km  and  the 
total  wind  at  RgaiesF  at  a  height  of  10  m  at  this  radius. 


An  alternative  depiction  of  the  vortex  evolution  is 
provided  by  the  Hovmoller  diagram  of  the  azimuthally- 
averaged  tangential  wind  component  near  the  top  of  the 
boundary  layer  (Fig.  5),  taken  nominally  at  a  height  of  1  km 
with  the  radius  of  the  wind  maximum  and  radius  of  gales 
at  this  level  indicated.  To  suppress  small-scale  noise  and 
highlight  the  main  features,  the  field  has  been  smoothed  by 
computing  a  3  h  time  average  centred  on  times  12  h  apart. 


(a)  time  (d) 


(b)  time  (d) 

Figure  3.  Time  series  of  azimuthally-averaged  quantities  at  latitudes 
10°N,  20°N,  and  30°N:  (a)  azimuthally-averaged  maximum  vertical 
velocity;  (b)  azimuthally-averaged,  12  hour  time  averaged,  radius  at 
which  the  maximum  vertical  velocity  at  a  height  of  6  km  occurs. 

This  figure  shows  also  the  progressive  expansion  of  the 
inner  and  outer  core  as  the  vortex  matures  and  decays  and 
will  prove  useful  in  the  discussion  of  the  boundary  layer 
behaviour  in  Section  4. 1 . 

3.3  Relationship  to  Smith  et  al.  (2011) 

The  foregoing  results  have  some  similarities  to  those  from 
the  axisymmetric  calculations  in  Smith  et  al.  (2011),  which 
were  designed  to  investigate  the  rotational  constraint  on 
the  intensity  and  size  of  tropical  cyclones  using  a  minimal, 
three-layer,  axisymmetric  model.  For  example,  the  vortices 
in  the  three  simulations  shown  here  have  a  similar  maxi¬ 
mum  intensity  during  their  mature  phases  as  do  those  in 
this  range  of  latitudes  in  Smith  et  al..  Moreover,  during  the 
mature  phase  in  the  calculations  here,  the  storm  at  20°N  is 
larger  than  those  at  10°N  and  30°N,  which  is  similar  also  to 
the  findings  of  Smith  et  al..  Flowever,  a  problem  with  such 
comparisons  is  that  the  vortices  in  the  Smith  et  al.  calcula¬ 
tions  reached  a  quasi-steady  state  after  a  few  days,  whereas 
those  in  the  present  study  ultimately  decay  in  terms  of  Vmax 
at  rates  that  are  dependent  on  latitude.  Thus  after  30  days, 
Vmax  at  10°N  is  appreciably  larger  than  that  at  30°N  and 
the  sizes  of  the  vortices  are  still  increasing.  As  shown  by 
Smith  et  al.  (2014),  the  existence  of  a  steady  state  requires 
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(a)  time  (d) 


0  3  6  9  12  15  18  21  24  27  30 


(b)  time  (d) 

Figure  4.  Time  series  of  azimuthally-averaged  radially  integrated 
mass  flux  at  latitudes  10°N,  20°N,  and  30°N  at  heights  of  (a)  1.5  km 
and  (b)  6  km.  This  quantity  is  calculated  by  integrating  the  vertical 
mass  flux,  where  it  is  positive,  radially  out  to  200  km,  and  then 
dividing  by  1  x  109  for  plotting  purposes. 

the  presence  of  a  steady  angular  momentum  source  in  the 
model  that  will  not  be  present  in  real  storms. 

3.4  Secondary  circulation 

The  strength  of  the  secondary  circulation  depends,  inter 
alia,  on  the  strength  and  radial  location  of  the  convective 
forcing,  which  are  controlled  by  the  dynamics  of  the 
boundary  layer.  Measures  of  this  strength  and  location 
are  the  maximum  azimuthally-averaged  vertical  velocity, 
u>max,  and  the  radius,  Rwrnax,  at  which  it  occurs.  The 
long  term  behaviour  of  these  quantities  is  shown  in  Fig. 
3.  The  largest  wmax  occurs  at  the  lowest  latitude  during 
the  intensification  phase  as  shown  in  SKM,  although  during 
the  30  day  integration,  wmax  is  similar  at  all  latitudes.  At 
20°N  and  30°N  the  largest  values  of  wmax  occur  in  the 
mature  phase,  and  the  location  of  wmax,  Rwmax,  changes 
little  during  this  period.  Subsequently,  wmax  decreases  at 
all  three  latitudes  and  Rwmax  moves  radially  outwards  with 
time,  the  radial  increase  decreasing  with  increasing  latitude. 

A  complementary  measure  of  the  strength  of  the 
azimuthally-averaged  secondary  circulation  is  the  verti¬ 
cal  mass  flux  carried  by  the  circulation  in  the  eyewall 
updraught.  This  quantity  characterizes  the  ability  of  the 
eyewall  updraught  to  ventilate  the  mass  flux  expelled 


20  N  MM5  vtop 


Figure  5.  Azimuthally-averaged  and  temporally-smoothed 
Hovmoller  plot  of  tangential  velocity  at  a  height  of  1  km  near 
the  top  of  the  boundary  layer  from  MM5  output  at  20°N.  Contour 
interval:  contours  every  5  m  s'1  with  the  17  m  s-1  contour  coloured 
blue.  The  interval  from  17  m  s^1  to  30  m  s^1  is  shaded  light 
red,  that  from  30  m  s-1  to  50  m  s-1  by  a  slightly  darker  red 
shading,  and  values  of  50  m  s-1  and  above  are  shaded  dark  red 
with  black  contours.  The  white  curve  represents  the  radius  of  the 
time  and  azimuthally-averaged  maximum  tangential  velocity  at  a 
height  of  1  km,  while  the  outer  blue  contour  represents  the  time  and 
azimuthally-averaged  radius  of  gales  at  this  height. 

by  the  inner-core  boundary  layer.  If  not  all  of  the  mass 
ejected  from  the  boundary  layer  can  be  accepted  by  eyewall 
updraught,  the  residual  will  flow  radially  outwards,  leading 
to  spin  down  of  the  tangential  wind  above  the  boundary 
layer  as  a  result  of  absolute  angular  momentum  conserva¬ 
tion.  Figure  4  shows  a  time  series  of  this  quantity  at  heights 
of  1.5  km  (panel  a)  and  6  km  (panel  b),  calculated  by  inte¬ 
grating  the  vertical  mass  flux  where  it  is  positive  radially 
out  to  200  km.  The  largest  mass  flux  at  both  these  heights 
occurs  at  the  lowest  latitude  throughout  the  30  day  inte¬ 
gration.  Interestingly  the  mass  flux  at  10°N  continues  to 
increase  during  the  decay  phase.  Thus,  even  though  the  eye- 
wall  expands  and  wmax  decreases  with  time,  the  amount 
of  mass  being  carried  by  the  updraught  at  a  height  of  6 
km  increases  at  this  latitude.  The  dependence  of  these  time 
series  on  latitude  will  be  discussed  further  in  section  6. 

3.5  Summary 

In  summary,  the  model  cyclones  grow  progressively  in  size 
and  decay  in  intensity  after  reaching  their  mature  stage.  In 
terms  of  either  Rgaies  or  RgaiesF ,  the  vortex  size  increases 
with  time  throughout  the  calculation,  except  at  the  largest 
latitude,  where  it  tends  to  become  steady.  The  strongest 
vertical  motion  occurs  during  the  intensification  phase,  and 
weakens  with  time,  while  the  location  of  the  strongest 
vertical  motion  moves  radially  outward  with  time  at  all 
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latitudes  studied.  The  evolution  of  vertical  mass  flux  carried 
by  the  eyewall  updraught  depends  also  on  latitude.  We 
explore  the  reasons  for  these  findings  in  section  6. 

4  Physical  interpretations 

The  tight  coupling  between  the  flow  above  the  boundary 
layer  and  that  within  the  boundary  layer  makes  it  impossi¬ 
ble  in  general  to  present  simple  cause  and  effect  arguments 
to  explain  vortex  behaviour.  The  best  one  can  do  is  attempt 
to  articulate  the  individual  elements  of  the  coupling,  which 
might  be  described  as  a  set  of  coupled  mechanisms.  In  this 
spirit  we  begin  in  the  next  subsection  by  examining  the  role 
of  the  expanding  outer  wind  profile  on  the  flow  within  the 
boundary  layer,  isolating  what  we  refer  to  as  a  boundary- 
layer  control  mechanism.  Then  in  subsections  4.2  and  4.3 
we  articulate  the  series  of  mechanisms  that  are  involved  in 
providing  boundary  layer  feedback  to  the  interior  flow. 

4. 1  The  boundary  layer  control  mechanism 

To  explore  the  behaviour  of  the  outer-core  expansion  of  the 
vortices  described  in  section  3  and  exemplified  in  the  case 
of  20°N  by  the  Hovmdller  diagram  of  tangential  wind  field 
shown  in  Fig.  5,  we  begin  by  examining  the  response  of 
the  boundary  layer  to  such  an  expansion.  To  this  end  we 
employ  a  simple,  steady,  slab  boundary  layer  model  as  in 
SKM.  The  details  of  this  model  and  the  justification  for 
its  use  are  given  in  SKM.  The  use  of  this  simple  model 
is  a  key  to  breaking  into  the  chain  of  coupled  mechanisms 
referred  to  above.  This  is  because  the  steady  boundary  layer 
equations  on  which  the  model  is  based  are  parabolic  in  the 
radially-inward  direction.  Thus,  the  inflow  and  hence  the 
ascent  (or  descent)  at  the  top  of  the  boundary  layer  at  a 
given  radius  knows  only  about  the  tangential  wind  profile 
at  larger  radii:  these  flow  features  know  nothing  directly 
about  the  vertical  motion  in  the  MM5  model  calculation 
at  the  top  of  the  boundary  layer,  including  the  pattern  of 
ascent  into  the  eyewall  cloud  associated  with  convection 
under  the  eyewall.  In  contrast,  the  MM5  model  does  not 
solve  the  boundary  layer  equations  separately  and  it  does 
not  make  any  special  boundary-layer  approximation.  Thus, 
the  ability  of  the  slab  boundary-layer  model  to  produce  a 
radial  distribution  of  vertical  motion  close  to  that  in  the 
time-dependent  MM5  model  provides  a  useful  measure  of 
the  degree  of  boundary-layer  control  in  the  evolution  of  the 
vortex. 

For  completeness,  a  brief  summary  of  the  slab 
boundary-layer  model  is  given  in  section  3  of  SKM.  For 
simplicity  we  assume  that  the  boundary  layer  has  a  con¬ 
stant  depth  of  1000  m,  comparable  to  the  depth  of  the  layer 
of  strong  inflow  in  the  MM5  simulation.  We  focus  now 
on  a  comparison  with  the  MM5  simulation  for  20°N.  The 
slab  boundary  layer  calculations  are  performed  every  12  h 
using  the  smoothed  tangential  wind  profile  extracted  from 
the  MM5  calculation  shown  in  Fig.  5.  Hovmdller  diagrams 
of  the  radial  and  tangential  velocity  components,  Ub,  i>& 


respectively,  and  of  the  vertical  velocity,  wu,  at  the  top  of 
the  boundary  layer  for  these  solutions  are  compared  with 
the  corresponding  diagrams  from  the  MM5  simulation  in 
Fig.  6.  The  radial  and  tangential  wind  components  from  the 
MM5  calculation  are  averaged  over  the  lowest  1  km  depth, 
corresponding  to  an  average  over  the  depth  of  the  boundary 
layer,  to  provide  a  fair  comparison  with  the  slab  boundary 
layer  fields. 

Despite  the  fact  that  the  slab  boundary  layer  calcu¬ 
lations  break  down  at  some  inner  radius  where  the  radial 
velocity  tends  to  zero  and  the  vertical  velocity  becomes 
large,  they  capture  many  important  features  of  the  corre¬ 
sponding  depth  averaged  boundary  layer  fields  from  the 
MM5  simulations.  For  example,  they  capture  the  broaden¬ 
ing  of  the  vortex  core  with  time,  i.e.  the  increase  of  Rvmax 
and  eyewall  location,  characterized  by  the  location  of  max¬ 
imum  Wh-  They  capture  also  the  broadening  of  the  outer 
radial  and  tangential  wind  field.  However,  they  overesti¬ 
mate  the  radial  extent  of  the  subsidence  outside  the  eyewall 
(compare  panels  (e)  and  (f)  of  Fig.  6).  For  reasons  artic¬ 
ulated  in  the  first  paragraph  of  this  section,  these  results 
provide  strong  support  for  the  existence  of  a  dynamical 
control  by  the  boundary  layer  on  the  evolution  of  the  vor¬ 
tex.  There  is  also  a  thermodynamic  control  that  we  explore 
in  the  next  section. 

4.2  Role  of  boundary  layer  thermodynamics  and  convec¬ 
tive  forcing 

The  boundary  layer  is  important  not  only  in  determining 
the  location  of  the  eyewall  and  the  radial  profiles  of  vertical 
velocity  and  tangential  momentum  entering  the  eyewall:  it 
plays  an  important  role  in  determining  also  the  radial  dis¬ 
tribution  of  diabatic  heating  within  the  eyewall.  The  reason 
is  that  the  wind  field  in  the  boundary  layer  affects  the  radial 
distribution  of  the  surface  enthalpy  flux.  In  turn,  this  distri¬ 
bution  is  important  because,  in  the  context  of  axisymmet- 
ric  balance  dynamics,  it  is  well  known  that  diabatic  heat¬ 
ing  associated  with  moist  deep  convection  in  the  eyewall 
leads  to  inflow  in  the  lower  troposphere  below  the  level  of 
the  heating  maximum  (Eliassen  1951,  Willoughby,  1979). 
Moreover,  the  strength  of  the  inflow  is  related,  inter  alia, 
to  the  negative  radial  gradient  of  diabatic  heating  rate.  This 
inflow  draws  M -surfaces  inwards  above  the  boundary  layer 
where  M  is  approximately  conserved  and  is  a  central  fea¬ 
ture  of  the  conventional  spin-up  mechanism  (Montgomery 
and  Smith  (2014)). 

We  show  in  Fig.  7  Hovmdller  diagrams  of  the  time 
and  azimuthally-averaged  diabatic  heating  rate  in  the  MM5 
simulations  at  a  height  of  6  km  for  all  three  experiments. 
Superimposed  on  these  plots  is  the  radius  of  the  time 
and  azimuthally-averaged  maximum  vertical  velocity  at  a 
height  of  6  km.  The  distribution  of  the  diabatic  heating 
rate  and  its  radial  gradient  weaken  considerably  with  time 
at  all  latitudes.  In  particular,  in  all  three  experiments,  the 
location  of  the  maximum  diabatic  heating  rate  and  the 
radius  of  maximum  vertical  velocity  move  radially  out¬ 
wards  with  time.  The  outward  movement  of  the  maximum 
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diabatic  heating  rate  is  expected  from  the  boundary  layer 
behaviour  described  above.  We  explore  the  consequences 
of  the  changing  heating  rate  for  the  calculation  at  20°N 
in  the  next  subsection.  The  latitudinal  dependence  is  dis¬ 
cussed  in  section  6. 

4.3  Diabatically-forced  balanced  overturning  circulation 

Following  on  from  the  foregoing  description  of  the 
behaviour  of  the  diabatic  heating  rate,  the  question  now 
is:  can  the  expansion  of  the  outer  wind  field  in  the 
MM5  calculation  be  explained  by  the  temporal  changes 
in  the  distribution  of  the  diabatic  heating  rate?  In  order 
to  answer  this  question,  we  adopt  the  approach  of  Bui 
et  al.  (2009)  and  solve  the  Sawyer-Eliassen  equation  for 
the  secondary  circulation  forced  by  the  distribution  of 
azimuthally-averaged  diabatic  heating  rate  derived  from  the 
MM5  output.  The  coefficients  of  this  equation  are  deter¬ 
mined  by  the  azimuthally-averaged  tangential  wind  from 
the  MM5  output  and  the  corresponding  balanced  tempera¬ 
ture  field  determined  using  the  method  described  by  Smith 
(2006).  The  tangential  wind  field  and  diabatic  heating  rates 
are  time  averaged  over  a  three  hour  period  centred  on  inter¬ 
vals  of  12  h  from  the  initial  time.  In  this  section  we  examine 
the  results  for  the  20°N  simulation. 

The  SE-equation  can  be  solved  only  if  it  remains 
elliptic  at  every  grid  point.  If  this  is  not  the  case,  a 
regularization  procedure  must  be  carried  out  to  restore  the 
ellipticity  at  such  grid  points.  Here  we  follow  the  ad  hoc , 
but  physically  defensible  method  suggested  by  Moller  and 
Shapiro  (2002),  which  is  described  in  more  detail  by  Bui 
et  al.  (2009). 

The  SE-equations  are  solved3  with  the  diabatic  heating 
rate  at  20°N  (see  Fig.  7b),  and  the  left  panels  of  Fig.  8  show 
Hovmoller  plots  of  the  Sawyer-Eliassen  solutions  for  the 
radial  and  vertical  wind  fields  at  heights  of  2.5  km  and  6 
km,  respectively.  The  right  panels  show  the  corresponding 
fields  from  the  MM5  simulations.  It  is  seen  that  the  SE- 
solutions  broadly  capture  the  magnitude  and  location  of  the 
eyewall  updraught  at  a  height  of  6  km.  There  are  differences 
in  detail  in  that  these  solutions  overestimate  the  inflow 
at  a  height  of  2.5  km  (panels  (c)  and  (d))  as  well  as  the 
horizontal  extent  of  the  region  of  subsidence  (panels  (a)  and 
(b)).  These  differences  would  be  consistent  with  the  neglect 
of  frictional  effects  in  the  SE-calculation. 

Figure  9  shows  Hovmoller  diagrams  comparing  the 
tangential  wind  tendencies  at  a  height  of  1.5  km  in  the 
balance  calculation  (excluding  friction),  in  both  cases  eval¬ 
uated  as  the  sum  of  the  radial  flux  of  absolute  vorticity 
and  the  vertical  advective  tendency.  Shown  also  is  Rgaies 
(pink  curve)  and  the  Rwmax  (red  curve)  as  a  function  of 
time.  Broadly  speaking,  in  the  MM5  calculation,  the  ten¬ 
dency  is  largely  positive  both  inside  Rgaies  and  outside,  at 

Specifically  we  solve  Eq.  (14)  of  Bui  et  al.  (2009)  neglecting  both 
frictional  forcing  and  the  relatively  small  contributions  of  the  "eddy 
terms”  as  defined  therein  in  this  equation. 


least  to  the  radius  of  400  km  shown.  However,  the  positive 
tendency  is  interspersed  by  intervals  of  negative  tendency 
that  can  be  traced  to  bursts  of  outflow  at  this  level  associ¬ 
ated  with  a  temporal  mismatch  between  the  boundary  layer 
convergence  and  the  ventilation  of  this  converging  air  by 
deep  convection  (not  shown).  The  corresponding  tendency 
in  the  SE-calculation  is  entirely  positive  with  values  at  and 
beyond  Rgaie s  increasing  in  strength  with  time.  The  lack 
of  negative  values  is  presumably  attributed  to  the  omission 
of  friction  effects  in  this  calculation  and  the  inability  to 
capture  the  nonlinear  dynamics  of  the  corner  flow  region 
(Smith  and  Montgomery,  2010). 

4.4  Synthesis 

In  subsection  4.1  we  showed  that,  forced  by  the  radially- 
expanding  azimuthally-averaged  tangential  wind  held  in 
the  outer  core  from  the  MM5  calculation  at  the  top  of  the 
boundary  layer,  a  simple,  steady,  slab  boundary  layer  pre¬ 
dicts  the  radial  expansion  of  the  radial  and  tangential  veloc¬ 
ity  components  in  the  inner  core  boundary  layer  as  well 
as  the  vertical  velocity  at  the  top  of  the  boundary  layer. 
The  radial  expansion  of  the  vertical  velocity  in  the  inner 
core  accounts  for  an  expansion  of  the  distribution  of  dia¬ 
batic  heating  rate  in  the  free  vortex  (subsection  4.2).  In  turn, 
we  showed  in  subsection  4.3  that  the  balanced  response  of 
the  vortex  to  the  expanding  azimuthally-averaged  diabatic 
heating  rate  diagnosed  from  the  MM5  model  leads  to  an 
expansion  of  the  inner  core  tangential  wind  held  at  the  top 
of  the  boundary  layer.  Taken  together,  these  hndings  pro¬ 
vide  an  explanation  to  the  question  posed  in  the  title  of  the 
paper. 

This  study  has  presented  a  system-scale  (azimuthally 
averaged)  diagnostic  analysis  of  vortex  evolution.  One 
notable  question  arises:  how  would  the  asymetric  eddies  be 
expected  to  inhuence  the  new  view  developed  here?  The 
“eddies”  as  vortical  hot  towers  play  an  important  role,  for 
example,  by  comprising  the  bulk  of  the  azimuthal  mean 
heating  rate  during  spin  up.  However,  at  later  times  when 
an  annular  eyewall  emerges,  the  towers  are  highly  strained 
as  discussed  in  Persing  et  al.  (2013).  It  is  not  immediately 
a  priori  clear  what  role  these  eddies  play  during  the  later 
expansion  of  the  vortex.  This  is  the  topic  of  a  future  study. 

5  An  alternative  explanation  for  size  expansion  in  two 
notable  tropical  cyclones 

The  foregoing  results  concerning  the  outward  migra¬ 
tion  of  the  tropical  cyclone  eyewall  accompanied  by  the 
slowly  declining  maximum  tangential  winds  provide  new 
insight  into  observations  of  mature  tropical  cyclones.  One 
well-documented  storm  that  shared  characteristics  simi¬ 
lar  to  these  idealized  numerical  experiments  was  category 
five  Atlantic  Hurricane  Isabel  (2003).  This  hurricane  was 
observed  with  multiple  reconnaissance  aircraft  over  three 
consecutive  days.  A  summary  of  the  inner-core,  vortex- 
scale  kinematic,  dynamic  and  thermodynamic  structure  of 
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Figure  6.  The  left  panels  show  azimuthally-averaged  and  temporally-smoothed  Hovmoller  plots  of  radial,  tangential  and  vertical  velocities  at 
the  top  of  the  boundary  layer  (a  height  of  1  km)  from  the  slab  boundary  layer  model,  with  a  constant  depth  of  1000  m  for  the  experiment  at 
20°N.  The  right  panels  show  the  corresponding  variables  front  the  MM5  output.  Contour  interval:  top  panel  contours  every  5  m  s-1  with  the 
17  m  s-1  contour  coloured  black.  The  interval  from  17  m  s-1  to  30  m  s-1  is  shaded  light  red,  that  from  30  m  s-1  to  50  m  s-1  by  a  slightly 
darker  red  shading,  and  values  of  50  m  s_1  and  above  are  shaded  dark  red.  Middle  panels:  first  contour  1  m  s^1,  then  in  5  m  s-1  intervals. 
Light  blue  shading  from  -5  m  s_1  to  -10  m  s_1,  darker  shading  below  10  m  s-1.  Bottom  panels:  first  (red)  contour  is  2  cm  s-1,  with  light 
red  shading  up  to  10  cm  s-1.  Darker  shading  from  10  cm  s-1  to  100  cm  s-1.  The  darkest  (red)  shading  is  enclosed  by  a  100  cm  s-1  black 
contour.  Regions  of  downward  motion  are  shaded  blue  and  enclosed  by  a  2  cm  s-1  blue  contour.  Solid  (red)  contours  positive,  dashed  (blue) 

contours  negative. 
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Figure  7.  Hovmoller  plots  of  the  azimuthally-averaged  and  temporally-smoothed  diabatic  heating  rate  in  the  MM5  simulations  at  a  height  of 
6  km.  Contour  interval:  1  K  hr-1  to  10  K  hr”1  shaded  light  red.  10  K  hr-1  to  25  K  hr-1  shaded  by  a  slightly  darker  red.  while  25  K  hr-1 
and  above  is  shaded  dark  red  with  black  contours  every  25  K  hr-1.  Shaded  blue  regions  are  below  -1  K  hr”1.  Solid  (red)  contours  positive, 
dashed  (blue)  contours  negative.  The  white  curve  represents  the  radius  of  the  time  and  azimuthally-averaged  maximum  vertical  velocity  at  a 

height  of  6  km. 


Isabel  over  these  three  days  of  observations  was  provided 
by  Bell  and  Montgomery  (2008).  One  particularly  intrigu¬ 
ing  and  unexplained  feature  of  Hurricane  Isabel's  evolution 
over  this  period  was  the  progressive  lateral  expansion  of 
the  eye  region  of  the  vortex  (see  Fig.  2  of  Bell  and  Mont¬ 
gomery  2008).  Figure  4  from  that  study  shows  that  the 
expansion  of  the  eye  was  accompanied  by  an  expansion  of 
the  radius  of  maximum  tangential  wind  ( Rvmax )  and  eye- 
wall.  In  particular,  the  Rvmax  grew  from  approximately  25 
km  on  12  September  to  45  km  on  13  September.  The  Rvmax 
continued  to  expand  to  55  km  on  14  September.  During 
this  extended  observation  interval,  the  maximum  tangential 
wind  declined  slowly  from  approximately  80  m  s”1  to  74 
m  s_1,  while  still  maintaining  a  category  five  status. 

As  discussed  by  Bell  and  Montgomery  (their  sec¬ 
tion  3a),  one  possible  explanation  for  the  expansion  of  the 
Rvmax  was  an  eyewall  replacement  cycle  between  12  and 
13  September.  Radar  reflectivity  imagery  and  microwave 
imagery  (their  Fig.  2a)  does  suggest  that  an  outer  rainband 
began  to  encircle  the  primary  eyewall  late  on  12  September. 
However,  limited  flight-level  data  and  microwave  imagery 
between  the  intensive  observation  periods  preclude  a  defini¬ 
tive  determination  of  an  eyewall  replacement  event.  The 
inner-core  expansion  mechanism  discussed  in  the  foregoing 
sections  offers  a  plausible,  and  arguably  simpler,  explana¬ 
tion  of  Isabel’s  inner-core  size  expansion. 

Another  documented  case  of  a  storm  that  intensified 
rapidly  with  no  eyewall  contraction,  but  a  significant  inner- 
core  size  increase  was  Typhoon  Megi  (2010).  On  the 
basis  of  a  numerical  simulation  of  this  storm,  (Wang  and 
Wang,  2013)  proposed  that  “the  inner-core  size  increase 
was  primarily  related  to  the  binary  interaction  of  Megi 
with  a  large-scale  low-level  depression  in  which  Megi 
was  embedded.  They  noted  that  “the  shearing/merging  of 
the  large-scale  depression  with  Megi  and  the  subsequent 
axisymmetrization  led  to  the  strengthening  of  the  outer 
circulation  of  Megi.”  As  shown  herein,  this  strengthening 


of  the  outer  circulation  would  lead  to  a  boundary-layer 
response  capable  of  explaining  the  inner-core  size  increase. 

6  Dependence  on  latitude 

The  foregoing  explanation  for  the  progressive  radial  expan¬ 
sion  of  the  vortex  is  based  on  a  calculation  for  20°N.  Here 
we  examine  the  latitudinal  dependence  of  this  intrinsic 
behaviour.  The  results  of  SKM  show  that  for  a  moderately 
strong  initial  vortex  (Vmax  =  15  m  s”1),  the  boundary  layer 
dynamics  exert  a  significant  control  on  the  vortex  spin-up 
rate  at  a  particular  latitude.  We  see  no  reason  to  believe  that 
this  control  should  diminish  in  prominence  as  the  vortex 
matures  and  decays.  Another  important  factor  is,  of  course, 
the  strength  of  the  convective  forcing. 

Some  of  the  latitudinal  behaviour  in  the  present  cal¬ 
culations  was  shown  in  Figs.  (1)  -  (4)  and  described  in 
section  3.  In  particular,  it  was  shown  that  during  the  ini¬ 
tial  period  of  contraction,  as  measured  by  rvmax,  the  mini¬ 
mum  value  of  rvmax  decreases  with  latitude.  As  the  vortices 
mature  and  decay,  rvmax  fluctuates  in  value  but  the  running 
mean  increases  as  the  difference  between  different  latitudes 
declines. 

During  the  first  half  of  the  simulations,  in  terms  of 
either  Rgaies  or  RgaiesF,  the  vortex  size  increases  above 
the  boundary  layer  with  time  at  all  three  latitudes  (Fig. 
2(a), (b)).  The  vortex  at  10°N  expands  a  little  more  slowly 
than  those  at  20°N  and  30°N,  which  have  essentially 
the  same  size  evolution.  During  the  second  half  of  the 
simulations,  the  vortices  at  10°N  and  20°N  grow  in  size 
together  while  the  size  of  the  vortex  at  30°N  becomes 
approximately  stationary,  except  for  a  small  increase  near 
the  end  of  the  simulation.  The  size  values  at  a  given 
time  are  not  uniformly  ordered  with  latitude  on  account 
of  two  competing  effects.  In  the  case  of  Rgaies >  these 
effects  are  the  convectively-driven  inflow,  which,  because 
of  changes  in  the  distribution  of  diabatic  heating  in  the 
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Figure  8.  The  left  panels  show  azimuthally-averaged  and  temporally-smoothed  Hovmoller  plots  of  radial  and  vertical  velocities  at  a  height 
of  2.5  km  for  radial  velocity,  and  a  height  of  6  km  for  vertical  velocity,  from  the  balanced  SE  equation  for  the  experiment  at  20°N.  The 
right  panels  show  the  corresponding  variables  from  the  MM5  output.  Contour  interval,  top  panels:  first  (red)  contour  2  cm  s_1,  with  light 
red  shading  up  to  10  cm  s-1.  Darker  shading  from  10  cm  s-1  to  100  cm  s-1.  The  darkest  (red)  shading  is  enclosed  by  a  100  cm  s-1  black 
contour.  Regions  of  downward  motion  are  enclosed  by  a  2  cm  s^1  dashed  blue  contour,  and  blue  shading.  Bottom  panels:  first  contour  0.5 
m  s-1  and  then  in  intervals  of  1  m  s-1  from  1  m  s^1.  Light  shading  from  1  m  s-1  to  2  m  s_1,  darker  shading  above  2  m  s-1.  Solid  (red) 

contours  positive,  dashed  (blue)  contours  negative. 


lower  troposphere,  decreases  in  strength  and  radial  extent 
as  the  latitude  increases  (Fig.  10),  and  the  radial  gradient  of 
absolute  angular  momentum,  which  increases  with  latitude. 
It  is  not  possible  a  priori  to  anticipate  which  effect  will 
dominate:  one  must  do  the  calculation  to  determine  the 
outcome. 

The  time  evolution  of  RgaiesF  is  similar  to  that  of 
Ftgaies  (compare  Figs.  2(a)  and  (b)),  but  the  latitudinal  dif¬ 
ferences  are  smaller.  In  this  case,  the  definition  of  RgaiesF 
includes  both  the  tangential  and  radial  velocity  components 
and  understanding  its  behaviour  requires  consideration  of 
the  radial  and  tangential  momentum  equations  including 
friction.  Again,  one  has  to  do  the  calculation  to  determine 
the  outcome.  Our  calculations  show  that  these  competing 


effects  nearly  cancel  with  only  small  differences  in  RgaiesF 
between  the  different  latitudes. 

We  showed  in  Fig.  3  that  the  strongest  azimuthally- 
averaged  vertical  motion  occurs  during  the  intensification 
phase,  where  it  is  notably  larger  at  10°N  than  at  the  higher 
latitudes.  Thereafter,  the  peak  values  weaken  with  time 
and  there  appears  to  be  little  systematic  difference  between 
them  in  the  range  of  latitudes  examined.  In  contrast,  the 
area-averaged  vertical  mass  flux  carried  by  the  azimuthal 
mean  eyewall  updraught  exhibits  rapid  temporal  fluctua¬ 
tions  with  time,  but  the  running  mean  increases  monotoni- 
cally  with  time ,  both  at  1 .5  and  6  km  (Fig.  4).  At  a  particular 
time,  the  mass  flux  is  less  when  the  latitude  is  increased. 
Figure  1 1  shows  the  difference  in  a  running  time  average 
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Figure  9.  Hovmoller  diagrams  comparing  the  azimuthally-averaged  and  temporally-smoothed  tangential  wind  tendencies  at  a  height  of  1 .5  km 
in  the  MM5  model  against  the  Sawyer-Eliassen  balance  calculation  (excluding  friction).  The  figure  shows  also  the  time  series  of  azimuthally- 
averaged  radius  of  gales  at  this  level  (pink  curve)  and  the  radius  of  maximum  vertical  winds  (red  curve).  Contour  interval:  0.5  x  10~°  m  s-2 
to  2  xlO-5  m  s'2  in  light  red  or  blue  shading,  dark  red  or  blue  shading  above  2  xlO-5  m  s-2  and  contours  in  intervals  of  2  xlO-5  m 
s~2  from  2  x  10-5  m  s~2.  Positive  tendencies  indicated  by  solid  red  contours  and  red  shading,  negative  tendencies  indicated  by  dashed  blue 

contours  and  blue  shading. 


30  N  MM3  u,  a  am  3.5  km 


Figure  10.  Azimuthally-averaged  and  temporally-smoothed  plots  of  AAM  and  radial  winds  at  a  height  of  2.5  km  for  (a)  ION,  (b)  20N  and  (c) 
30N.  Contour  intervals:  AAM  contoured  in  black  from  2  m2  s-1  in  intervals  of  2  m2  s_1.  Radial  velocity  contours  at  0.2  m  s_1  and  0.5  m 
s_1.  Light  shading  between  0.2  m  s^1  and  0.5  m  s-1,  darker  shading  above  0.5  m  s-1  with  the  darkest  shading  above  1  m  s_1.  Solid  (red) 
contours  positive,  dashed  (blue)  contours  negative,  red  shading  positive,  blue  shading  negative. 


over  12  hours  of  the  corresponding  mass  flux  curves  in  Fig. 
4  between  6  km  and  1.5  km,  expressed  as  a  percentage  of 
the  mass  flux  at  1.5  km.  During  the  rapid  intensification 
phase  (approx.  1-6  days  depending  on  latitude),  when  the 
convection  is  most  vigorous  (Fig.  3  a),  there  is  adequate 
ventilation  of  the  air  that  exits  the  boundary  layer  at  all  lati¬ 
tudes,  in  the  sense  that  the  mass  flux  at  6  km  height  exceeds 
that  at  1.5  km. 

As  the  vortices  mature  (approx,  during  the  period  3- 
9  days  depending  on  latitude),  the  boundary  layer  upflow 
slightly  exceeds  the  mass  flux  at  6  km,  again  at  all  latitudes. 
This  behaviour  is  accompanied  by  a  progressive  warming 
of  the  upper  levels  (not  shown),  which  reduces  the  convec¬ 
tive  instability  in  the  inner  core  region  as  is  evident  in  the 


progressive  decline  in  maximum  azimuthally-averaged  ver¬ 
tical  velocity.  The  increase  in  the  mass  flux  at  6  km  beyond 
the  intensification  phase,  despite  the  reduced  convective 
instability,  is  accompanied  by  an  outward  radial  displace¬ 
ment  of  the  azimuthal  mean  eyewall  updraught  as  well  as 
a  broadening  of  the  updraught.  During  the  decay  phase 
(beyond  about  9  days),  the  mass  flux  at  a  height  at  6  km 
once  again  exceeds  that  at  1 .5  km  and  is  accompanied  by  a 
continued  structural  change  in  the  eyewall. 

The  structural  evolution  of  the  eyewall  updraught 
referred  to  above  is  illustrated  in  Fig.  12,  which  shows 
verical  cross-sections  of  the  azimuthally-averaged  vertical 
velocity  in  the  calculations  at  10°N  and  30°N  at  selected 
times.  Shown  also  is  the  5ms-1  contour  of  radial  velocity 
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time  (d) 

Figure  11.  Normalized  mass  flux  difference  between  6  km  and  1.5 
km  for  10-30N.  The  difference  is  calculated  as  a  running  time  average 
over  12  hours  of  the  corresponding  mass  flux  curves  in  Fig.  4  between 
6  km  and  1.5  km,  and  then  expressed  as  a  percentage  of  the  mass  flux 
at  1.5  km. 

and  the  slope  of  the  eyewall  in  the  lower  troposphere, 
characterized  by  the  difference  in  radial  location  of  wmax 
between  heights  of  1 .5  km  and  6  km.  At  6  days,  when  the 
two  vortices  have  reached  maturity,  the  eyewall  updraught 
is  upright  at  both  latitudes,  but  in  comparing  the  details 
of  these  figures  one  should  recall  that  the  vortex  at  30°N 
reaches  maturity  three  days  later  than  that  at  10°N.  A 
prominent  difference  between  the  two  simulations  is  the 
more  inward  penetration  of  the  boundary  layer  inflow  at 
the  lower  latitude  and  the  correspondingly  larger  region  of 
outflow  just  above  the  inflow  layer.  This  more  extensive 
outflow  at  10°N  appears  to  influence  the  radial  location  and 
extent  of  the  updraught  at  all  levels. 

As  time  proceeds,  the  updraughts  at  both  latitudes 
continue  to  broaden  as  the  pattern  of  upflow  at  the  top  of 
the  boundary  layer  broadens  (Fig.  13)  and  the  maximum 
vertical  velocity  aloft  weakens  (Fig.  12).  This  broadening 
is  accompanied  by  a  progressively  increasing  outward  tilt 
of  the  updraught  as  indicated  in  Fig.  12.  Moreover,  the 
positive  ventilation  implied  by  the  larger  mass  flux  at  a 
height  of  6  km  than  at  1.5  km  is  accompanied  by  (mostly) 
entainment  into  the  updraught  at  outer  radii  as  seen  in 
the  Flovmoller  diagram  of  inflow  and  the  radially-inward 
slope  of  the  M-surfaces  a  height  of  2.5  km  in  Fig.  10.  A 
consequence  of  this  inflow  is  a  continued  expansion  of  the 
outer  tangential  wind  field  in  the  lower  troposphere.  Note 
that  at  30°N,  the  M-surfaces  at  2.5  km  height  beyond  a 
radius  of  about  300  km  and  after  about  15  days  are  no 
longer  moving  inwards,  indicating  that  the  outer  wind  field 
has  ceased  to  expand,  consistent  with  the  behaviour  of 
Rgaies  at  this  latitude  in  Fig.  2a. 

In  summary,  the  purpose  of  this  section  was  to  describe 
and  interpret  the  variability  of  tropical  cyclone  lifecycle 
with  latitude  in  a  quiescent  environment.  We  have  shown 
that,  at  all  latitudes,  the  tangential  wind  field  in  the  free 
troposphere  progressively  broadens  on  account  of  inflow 
associated  with  deep  convection  in  the  eyewall  updraught, 
except  at  30°N,  where  the  broadening  declines  after  about 


15  days.  The  inflow  is  consistent  through  mass  continuity 
with  the  excess  mass  flux  in  the  eyewall  updraught  at  a 
height  of  6  km  compared  with  that  at  a  height  of  1.5  km 
during  the  decay  stage  of  the  vortices.  The  convective  mass 
fluxes,  both  at  6  km  and  1.5  km,  increase  with  decreasing 
latitude,  a  feature  that  may  be  attributed  to  the  strength 
of  the  boundary  layer  control  on  the  eyewall  updraught, 
which  increases  also  with  decreasing  latitude.  In  brief,  the 
boundary  layer  affects  the  location  and  radial  extent  of 
the  eyewall  updraught  as  well  as  the  radial  gradient  of  de 
exiting  the  boundary  layer  into  the  updraught.  As  shown 
above  and  in  SKM,  these  features  of  the  boundary  layer 
depend  both  on  the  radial  distribution  of  tangential  wind  at 
the  top  of  the  boundary  layer  as  well  as  on  latitude.  Because 
of  the  tight  coupling  between  the  boundary  layer  and  the 
tangential  wind  field  at  its  top,  it  is  not  possible,  apriori,  to 
make  a  prediction  about  which  effect  will  win:  one  has  to 
do  the  calculation. 

7  Conclusions 

We  have  examined  the  long  term  (out  to  30  d)  behaviour 
of  tropical  cyclones  in  the  prototype  problem  for  cyclone 
intensification  on  an  /-plane  using  a  nonhydrostatic  three- 
dimensional  numerical  model.  We  showed  that  the  model 
cyclones  grow  progressively  in  size  and  decay  in  intensity 
after  reaching  their  mature  stage.  As  they  decay  the  inner 
core  size,  characterized  by  the  radius  of  the  eyewall,  and 
size  of  the  outer  circulation,  characterized  for  example  by 
the  radius  of  gale  force  winds,  both  progressively  increase, 
except  at  latitude  30°N  where  the  radius  of  gales  tends  to 
become  steady.  The  strongest  vertical  motion  occurs  during 
the  intensification  phase,  and  weakens  with  time,  while  the 
location  of  the  strongest  vertical  motion  moves  radially 
outward  with  time  at  all  latitudes  studied.  The  evolution  of 
vertical  mass  flux  carried  by  the  eyewall  updraught  depends 
also  on  latitude. 

We  answer  the  question  posed  in  the  title  of  the  paper 
by  isolating  what  might  be  described  as  a  set  of  coupled 
mechanisms,  the  first  of  which  we  refer  to  as  a  boundary- 
layer  control  mechanism.  We  showed  that,  forced  by  the 
radially-expanding  azimuthally-averaged  tangential  wind 
field  in  the  outer  core  from  the  numerical  model  simula¬ 
tion  at  the  top  of  the  boundary  layer,  a  simple,  steady,  slab 
boundary  layer  model  correctly  predicts  the  radial  expan¬ 
sion  of  the  radial  and  tangential  velocity  components  in  the 
inner  core  boundary  layer  as  well  as  the  vertical  velocity  at 
the  top  of  the  boundary  layer.  The  radial  expansion  of  the 
vertical  velocity  in  the  inner  core  accounts  for  an  expan¬ 
sion  of  the  distribution  of  diabatic  heating  rate  in  the  free 
vortex.  In  turn,  we  showed  that  the  balanced  response  of 
the  vortex  to  the  expanding  azimuthally-averaged  diabatic 
heating  rate  diagnosed  from  the  numerical  model  leads  to 
an  expansion  of  the  inner  core  tangential  wind  field  at  the 
top  of  the  boundary  layer.  Taken  together,  these  findings 
provide  an  explanation  to  the  question  posed  in  the  title 
of  the  paper.  They  provide  also  a  plausible,  and  arguably 
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Figure  12.  vertical  cross  sections  of  the  azimuthally-averaged,  24  hour  time  averaged  vertical  velocity  at  10N  and  30N.  Contour  interval: 
contours  0.1  m  s'1  to  0.9  m  s_1  enclosed  by  light  shading.  Thick  contour  above  1  m  s-1  encloses  an  area  of  darker  shading.  Regions  of 
downward  motion  are  shaded  blue  and  upward  motion  shaded  red.  The  black  contours  show  regions  of  strong  radial  inflow  and  outflow, 
greater  than  5  m  s-1.  Solid  contours  positive,  dashed  contours  negative.  Shown  also  is  the  slope  of  the  eyewall  in  the  lower  troposphere  (blue 
line),  characterized  by  the  difference  in  radial  location  of  WmOX  between  heights  of  1.5  km  and  6  km. 


Figure  13.  Hovmoller  plots  of  the  vertical  velocity  in  the  MM5  simulations  at  a  height  of  1.5  km.  Contour  interval:  first  (red)  contour  is  2  cm 
s_1,  with  light  red  shading  up  to  10  cm  s'1.  Darker  shading  front  10  cm  s_1  to  100  cm  s_1.  The  darkest  (red)  shading  is  enclosed  by  a  100 
cm  s_1  black  contour.  Regions  of  downward  motion  are  shaded  blue  and  enclosed  by  a  2  cm  s_1  blue  contour.  Solid  (red)  contours  positive, 

dashed  (blue)  contours  negative. 


simpler,  explanation  for  the  expansion  of  the  inner-core  of 
Hurricane  Isabel  (2003)  and  Typhoon  Megi  (2010)  during 
and  after  peak  intensity  than  given  previously. 

Finally  we  showed  that  the  strength  of  the  boundary 
layer  control  on  the  location  and  radial  extent  of  the  eyewall 
updraught,  as  well  as  the  thermodynamic  control  of  the 
boundary  layer  on  the  distribution  of  diabatic  heating  in 


the  updraught,  both  have  a  latitudinal  dependence.  Since 
the  distribution  of  diabatic  heating  determines  primarily 
the  strength  of  the  secondary  circulation,  it  has  an  indirect 
effect  on  the  radial  distribution  of  tangential  wind  at  the 
top  of  the  boundary  layer  which  therefore  has  a  latitudinal 
dependence  as  well.  Because  of  the  tight  coupling  between 
the  boundary  layer  and  the  tangential  wind  field  above  the 
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boundary  layer,  which  is  itself  determined  by  the  diabatic 
heating,  it  is  not  possible  to  make  an  a  priori  prediction  of 
how  the  vortex  evolution  will  change  with  latitude:  one  has 
to  do  the  calculation. 
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